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ABSTRACT. In this study we display a process for inverting the generalized Van- 
dermonde matrix, using the analytic properties of a fundamental system related 
to a specific linear difference equations. We establish two approaches allowing 
us to provide explicit formulas for the entries of the inverse of the generalized 
Vandermonde matrices. To enhance the effectiveness of our the approaches, sig- 
nificant examples and special cases are given. 
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1. INTRODUCTION 


The usual Vandermonde systems of equations of order r is given by, 
F 

SS NES 05. R= ease ~ 1 (1) 

i=1 
where the x; (1 < i < r) are the unknown variables, the A; (1 < i < r) are distinct real 
(or complex) numbers and the vn (0 < n < r — 1) given real (or complex) numbers. Let 
m; >1(1<1i<-s) bes integers and \; (1 < i < s) be distinct real (or complex) numbers. 
For a given real (or complex) numbers v, (0 < n < r — 1), where r = mı + --- + Ms, the 
related generalized Vandermonde systems of equations is defined as follows, 


s m,—-1 
J ) Tijn | AF SH thes n=0,1,...,r— 1, (2) 
i=1 \ j=0 


where the xj; (1 < i < s, 0 < j < m; — 1) are the unknown variables. The generalized 
Vandermonde system of equations is also known in the literature as a nonsingular usual 
Vandermonde system. The generalized Vandermonde systems (1)-(2) appear in several 
topics of mathematics such that the linear algebra, numerical analysis and polynomial 
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approximation or interpolation. They also own several important applications in var- 
ious areas of applied sciences and engineering like signal processing, statistics, coding 
theory and control theory (see for example [8, 10,14], and references therein). In order 
to solve Vandermonde systems (1)-(2), several methods have been provided in the lit- 
erature (see, for instance, [2, 6, 11, 12,16, 20]). It was shown that solving these systems 
is related to the inverses of their associated matrices, called also usual and generalized 
Vandermonde matrix, respectively. For inverting the usual Vandermonde matrix, several 
methods have been considered in various studies (see [2, 6,9, 11,17,19], and references 
therein). The search for an efficient approach for computing the inverse of the general- 
ized Vandermonde matrix, is still an attractive topic, because of its interest in various 
topics of mathematics and applied sciences. Especially, two recent methods have been 
elaborated in [2,12], for solving the generalized system (2). The process of [2] is base 
on the analytic formula (or the so-called analytic Binet formula) of the linear recursive 
sequences {vn }n>0 of constant coefficients and order r, defined as follows, 


Un = G0Un—1 + Q1Un—2 + +++ + ar—1Vn-r for every n > r (3) 
Un = An for n=0,1,...,r—1, (4) 


where apo, ..., ar—ı are the coefficients and ag, ..., &r—1 are the initial data. That is, the 
analytic formula of sequence (3)-(4), is given by, 


s mi—1 
Un = > bA Bignt | A2, (5) 
i=1 \ j=0 


where the A; (1 < i < s) are the roots of the (characteristic) polynomial P(z) = z" — 
agz”—!—--+-—ap_1, of multiplicities m1, m2,..., Ms, respectively. The scalars b; j (1 < i < s, 
0 < j < m; — 1) are obtained by solving a linear system, 


s mi—l1 
`> Toan A =n, O<n<r-l. (6) 
i=1 \ j=0 


We show that the system of equation (6) is nothing else but a generalized Vandermonde 
system (2). Indeed, the unknown variables x; j (1 < i < s,0 < j < m; — 1) are identified 
to the scalars §;,; (1 < i < s, 0 < j < m; — 1), namely, the two systems (2) and (6) are 
identical. Therefore, exhibiting explicit compact analytic formulas (5) for the general term 
Un, namely, explicit formulas for the scalars f; j (1 < i < s,0 < j < m;—1) without solving 
the generalized Vandermonde system (2), will permit us to establish explicit formula for 
the solution of this system, or equivalently, the entries of the inverse of the associated 
generalized Vandermonde matrix. 


This study concerns the implantation of a process for computing explicit formulas for 
the entries of the inverse of generalized Vandermonde matrix, through knowledge of the 
analytic formulation (5) of the fundamental system related to Expression (3), considered 
as a difference equation. To reach our goal, we exploit the recent studies of [1,3] in order 
to set up two approaches which make it possible to highlight a new explicit form of the 
analytical formula (5) of vn, namely, we are going to exhibit new explicit formulas for 


the scalars 6; (1 < i < s, 0 < j < m; — 1), and by the way we set out the inverse 
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of the generalized Vandermonde matrix, or equivalently we solve the the generalized 
Vandermonde system (2). In order to better understand our results, we will exhibit some 
special cases and illustrative examples. Finally, we analyze and discuss the results issued 
from our two approaches, by comparing them with the literature, especially with results 
of [2,12]. 


The remainder of this paper is organized as follows. Section 2 in devoted to the prop- 
erties of the fundamental system of Equation (3) and its related dynamic solution. Some 
explicit compact analytical formulas of the dynamic solution are proposed. Section 3 con- 
cerns the process of construction of the inverse of the generalized Vandermonde matrix, 
with the aid of the analytical expression of the fundamental system of Equation (3). The 
two main results on the inverse the generalized Vandermonde matrix related to the the 
generalized Vandermonde system (2), are provided in Section 4. Analysis and discussion 
are considered in section 5. Finally, conclusion and perspectives are given in Section 6. 


2. EXPLICIT ANALYTIC SOLUTION OF THE DYNAMIC SOLUTION OF (3) 


2.1. Fundamental system and its dynamic solution. Let Ef) (ao, ..., @r—1) be the K-vector 
space of finite dimension r, of solutions of Equation (3) of coefficients ag, ..., @r_1 Æ 0. 
Let HoP } nso; 0 < p < r — 1} be the family of sequences of EL (ao, ..+;@p—1), indexed 
by p (0 < p < r — 1) defined as follows, 


wP = aov? F av), eee Orv”), for n >r, (7) 
yl?) = p,n forO<n<r-l. 


Namely, this family {fu 0! 0 < p < r — 1} represents r copies of sequences (3) with 

mutually different sets of initial values, viz. uP = pn (O0 <n<r-1,0<p<r-!1), 

where 6,7, is the Kronecker symbol. For every {vn}n>o in EL (ao, ..+,@,—1) Of initial 
r—1 

data ao,..., @r—1. Let {Wn}n>o be the sequence defined by wn = `> a, uv), Then, as 
p=0 

a linear combination of the toth nso, we can show that {w,}n>0 is in E (aa, see ey Ny 


moreover we have wk = ax, for 0 < k < r—1. Therefore, using Lemma 2.1, we can 
r-1 


verify that vn = wn, for every n > 0, namely, vn = X apy’, for every n > 0. On 


p=0 
r-1 


the other side, suppose that >D a,v?) = 0, for every n > 0. Since uP) = Onn for 0 < 


p=0 
r—1 


n < r— 1, we show that `> apu) =a, = 0 for 0 < p < r — 1. Therefore, the family 
p=0 

HoP } nso; 0 < p < r — 1} is a basis of the K-vector space EM (ao, ...,@p—1). The family 

{{v? 103 0 < p < r—1}is known in the literature as the fundamental system of Equation 

(3), or the Fibonacci fundamental system of Equation (3). 


For reasons of utility, we first state the following elementary lemma, concerning the 
equality of two elements of the space ef (ao, ...,@r—1). 


Lemma 2.1. Let {Un}n>0 and {Wn}n>0 be in Ef (ao, eiri „ar—1). If Uk = Wk, for 
k =0,1,...,r — 1, then, we have vn = wn, for every n > 0. 


The proof of this lemma is done by simple reasoning by induction. 


Among the element of the fundamental system (7), the next theorem shows the closed 


relation between {oY 1.50 and the other elements {uP MY for 
p=0,...,7r—2. 


Theorem 2.2. Let HoP } nso; 0 < p < r—1} be the fundamental system of Equation (3). Then, 
for every p, with p = 0,...,r — 2, we have, 


vP) = Opies ies pur eu te + Ap yo ag (8) 


for every n> r. 


Proof. Let first recall that the sequences {ul MS are defined by (7), namely, 


oo .r—1, 
v 


DE cat® 


go eee ar10, forn >r. 


Let {wo} nso be the sequence defined by wi ) = and w = a, T D. We show that 


{wo} aso satisfies Equation (3), with initial data, 


(r— n 


w = 1 and w® = ap 1u,_1 =0, foreveryl<n<r-1. 
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Therefore, using Lemma 2.1, we derive wp’ = = th for every n > 0, namely, we have 


yO) = ar10 i D, for every n > 1. For p = 1, let wP }n>0 be the sequence defined by 


(1) (r— a) (r—1) 


we’ =9, wo) = 1 and w® =a, 2U 1 AUS 
It is clear that wP = diata J w®,, thus the sequence {uf DAN satisfies Equation 
(3), with initial data, 


w =0, w = 1 and w® =a,_ vt P wO, =0fr2<n<r-1. 


Therefore, Lemma 2.1 shows that the two sequences uP n>0 and {y) }n>0 are identi- 


cal, namely, we have aUe = Ap— ou Pe, + Gp_1U ol A More ans with the aid of the 
similar argument, for 2 < p < r — 2, we ee R sequence {wt ARR defined by the 
initial conditions wP ER wt) i=0, wP l= 1, and 

wP) = E E Ydre a ey Bape sy ww eee 
We can observe that wi?) = a, ee hi D4 wt- ae , where we Su = Opp "5? +e + 
Ap— or Da oM Hence, the sequence {w?) es ee Equation (3), with initial data, 
wP) mae we), 0, wP) = 1 and w®P) = a,_ eae yo tw? = 0 for p+1 < n < r-1. 


Hence, by applying Lemma 2.1, we derive that the two sequences {u }n>0 and P }n>0 


are identical. Therefore, for every p (0 < p < r — 2), we have 
vP) = app- TNE nA pur ae te + Gp wo a 
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The result of Theorem 2.2 can be also established by induction on p. However, we have 


used here an elementary process based on Lemma 2.1 and the fact that Ef) (a0, ...,ar—1) 
is a K-vector space 


Expression (8) shows that the sequence CTT will play a central role in the se- 
quel. The sequence a i is called dynamic solution of Equation (3). 


For illustrative purpose, we propose the following special case. 
Example 2.3. For r = 4, Expression (3) takes the form 


Un = a0Un—1 + A1Un—2 + a2Un—3 + a3Un—4 for n> 4. 
Therefore, the terms yO), yD), vP and v of the fundamental system are expressed in terms of 
. p (r-1) _ (3) 
the dynamic solution vp ” = vn 


G, vy) = azv, st azv, y?) = ayv®), F azv, oe azv®) g, 


under the form, 
of = a3U 
forn > 4. 


2.2. Analytical formulas of the dynamic solution: First approach. Let à; (1 < i < s) 
be the roots of the characteristic polynomial P(z) = z" — agz’~! — --- — ap—ı, related 
to sequence (3), of multiplicities m1, m2,..., Ms, respectively. For every m; > 1 we set 
ell = {(n1, byg ns) E€ NsTt: ni +e + Ni- + Nnipi +: Nns = mi — k — Tk In [5, Section 
4.1], the following expression was considered, 


G + M3 — ‘) 
Honos T A o 
p (A1; -s Às 


.— ).)\nj+tm; 
ell \ 1siAi<s Ar) 


for0<k<m,—land1<i<-s. 


Example 2.4. Let r = 7 and suppose that s = 3, mı = 2, ma = land mg = 4. Fori=3,k=1 


we have ra = {(n1,n2); nı + ng = 4 — 1 — 1 = 2}. Therefore, we have, 


Ce (aa) 
3] ~(-1)3 ie z 
Yi (A1, A2, A3) = (—1) (Ài = A3)mtma f (AQ = Ag)r2atme ` 


nı +n2=2 


nı +1 
nı 


Therefore, since és = land ( 
n2 


) = nı + 1l, we obtain, 


nı +1 1 
(Ài = A3)™+2 (2 = Ag)r2ti É 


af (Ar, Az, As) = (1) 
nytng=2 
The analytic formula of the dynamic solution ul’) related to the fundamental system 
(7), can be expressed in terms of the roots A; (1 < i < s) of the polynomial P(z) = 
zT mya A Bee ay by the previous expressions (9) of the On, ...,;As). More 
precisely, we have the following result. 
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Proposition 2.5. The analytic of the dynamic solution of) is given by the following 


ml 
formula v? ay (© y oe Pat) dj, for every n > r, where 


4=1, k=0 
ty Sy HOAs) 
Pik =g s(t, a ne (10) 
t=k i 


and the s(t, k) are the Stirling numbers of the first kind. 


s m,—-1 
Indeed, from [3, Theorem 2.2], we have yl Y =e - 6 HO, .. >) me 


i=1 
for all n > r. On the other hand, it is well known that 


C = n(n— 1)---(n— k+ 


k 
[j= `> s(k,t)nt, where the s(k,t) are Stirling numbers of the first kind. Therefore, the 
t=0 


combinatorial expression o) can be also expressed in terms of the Stirling numbers of 


k 
A : n\ _ areal e- (n-k+1) s(k,t) 4 
the first kind as follows ( À = il = 2 p Thus, we get 


s mi—l k [i] À À 
Tae > ( F (>: s(k, on) rea) Aç, for every n > r, or equivalently, 


i=1 k=0 t=0 


s mi—1 mi—l [i] 
yD) = 5y ( 5y ( z s(t, yey ee) n) Az, for every n >r. 


i=1 k=0 t=k 


Therefore, the results follows, namely, Expression (10) is established. 


More generally, the result of Proposition 2.5 allows us to determine the expressions of 
vE- a (1 < d< r -— 1). In summary, we have the following result. 


Proposition 2.6. Under the preceding data, ford = 1,...,r — 1, we have 


s mi—l1 
wp = oS 5 ot) A; for every n >r, 


i=1 \ j=0 
where 
mi—1 
i eee k 
JEE E (0) a, 1 
Cs 2 (D Pik j ( ) 


such that the foe are as in (10). 
Indeed, from Proposition 2.5, we derive, 
s mi—1 mi— l1 i 
wa =D ( J Ben - a) etsy. 2 Bie” > (nit ard |e. 
j=0 


i=1 k=0 i=1 
111 


s mi—l1 mi—1 

Therefore, we get v of ) = 5 ( ($ > Os uel he) e) Xi, which 

i=1 \ j=0 

allows us to derive Expression (11). 
For reason of clarity, in the following corollary, we consider the special useful case, for 


illustrating Propositions 2.5 and 2.6. 


Corollary 2.7. Special case s = 2: mı > 2 and mo > 2. Under the data of Propositions 2.5 
and 2.6, suppose that s = 2: mı > 2 and mo > 2, then, we have, 


mı—1 m2—-1 
-1) _ & acta) af + (35 ae! 3 DE n 


for every n > r = mı + Mg, where the BM Sa i) i = 1, 2areas in (10) , namely, 
mı—1 [1] m2-1 [2] 
Pa ye (A1, À ee y (AL, A 
(i. SL eG ye an and pS = So s(t, yw (12) 
t=k l t=k i 


and the s(t, k) are the Stirling numbers of the first kind. Moreover, for every d = 1,...,r — 1, we 
have, 


for every n > r = mi + mo, where the ee 


mı—l1 
o =A? `> (—1) ys oes ae Je j and C$ ) ae d osc 1)*- Pe C) dEi, 
= 


k=j 


are as in (11), namely, 


such that the be D and BE D are given by (12). 


2.3. Analytical formulas of the dynamic solution: Second approach. For a given se- 
quence (3), it was shown in [13] that, 


Un = p(n, r) Ao + p(n —1,r)Ar+-+--+ p(n -=r + 1,r)Arı, (13) 
for every n > r, where Aj = ar—1vi +--+ + Givp_1 (0 < i < r—1) and 
(ko +ki +-+ ee Lae ki kra 


p(n,r) A ko! ky! +++ kp—1! ao ay Qr—1 ’ ( ) 
ot2kit+--+rkp-1=n—-r 


with p(r,r) = 1 and p(n,r) = 0 for n < r—1. Let à; (1 < i < s) be the roots of the polyno- 
mial P(z) = 2” — aoz"! — - - - — ap—1, of multiplicities m1, m2,..., Ms, respectively. Using 
the divided difference techniques and Newton interpolation method, it was established 
in [1, Theorem 3.1] that Expression (14) of p(n, r), can be formulated in terms of the roots 
Ai (1 < i < s) and their multiplicities m1, m2,..., Ms as follows, 


fin Os) 
pir) = bp are for every n > r, (15) 
i=1 
112 


n-1 


with p(r,r) =1, p(n,r) =0 for O< n<r-—1,and fin(x) = — z , where 
I] (@-Ak)™* 
k=1, k#i 
ie (x) means the derivative of order k of the function fin. Especially, when the roots A; 
(1 < i < s) are simple, namely, s = r and mı = mz = --- = Mmr = 1, Expression (15) takes 
the form, 
s pD (\;) r ane 
p(n,r) = =D = ‘> 3 , (16) 
i=1 : © #10 J (A — Ak) 
k=1, k#i 


Expression (16) has been also established in [2,4]. In addition, let {wy }n>o0 be the se- 
quence defined by wn = p(n + 1, r), then Expression (13) shows that {wn }n>0 satisfies the 
recursive relation (3), and its initial conditions are wo = --- = wr—2 = Oand w;_; = 1. 


Hence, the dynamic solution meee eer and the sequence {wn }n>0 Satisfy the same re- 
cursive relation (3) and own the identical initial conditions. Therefore, Lemma 2.1 shows 
that, 


uv) = p(n + 1,r), for every n > 0. (17) 


n 


Therefore, taking into account Expressions (15) and (17), we derive that the analytic for- 
mula of the dynamic solution is given by, 


pe foo Os) 
ee Scr 


i=1 


for every n > r, (18) 


where the function fin+1 (1 < i < s) are defined as by fin41 (x) = 


k= 


= 
>” 
KK 
>. 


For m; = 1 we have m;—1 = 0, therefore, we get EPE (£) = fin+ı (£) = E 


For technical reasons, we set fin+1(£) = dn(x)Hi(x), where qn(£) = z” and H;(x) = 
1 


=a 


ne 


k=1, ki 


For m; > 2, let compute the explicit formula of the derivative of order m; — 1 of the 
g” 


function fin+1ı(£) = . To achieve our goal, we will proceed in two 


k=1, k#i 
steps. First, let f, g two functions admitting derivatives of order m > 1 on a nonempty 


subset of R. It is well known that, we have (f g)™ = >. o fO gd, Application of 
d=0 
this former formula to fin+1ı(x) = qn(x)H;(x), we obtain, 


aO = (PPE = (7) aaa), 


d=0 


where (n)a = n(n — 1) --- (n — d + 1). Since (n)a = ya s(d,h)n", where the s(d, h) are 
the Stirling numbers of the first kind, we show that, 


d=0 d=0 h=0 
m d m m 
Using the identity > > Tdh = X > Lad,n, for a bi-indexed sequence £a,n, we obtain 
d=0 h=0 h=0 d=h 


f(a) = D (È s(d,h) e sI ore) n'a”. Therefore, for m; > 2 we set £ = 


h=0 \d=h 
A; and m = n; — 1, thus we arrive to have, 


mi—l /mi—1 
Mi— - f a, 1 m,—d— = 
ier On) =e & s(d, an Ja AD JA ‘ nh. (19) 


h=0 d=h 


Second, to improve Expression (19), we will give the explicit form of the p — th derivative 
of the function H;(x). For this purpose, we use the following well-known formula, 


(fi for: fe) = iD e k Ta 
se ks) 35 


m m! . f i 
where @ “ | ie ail and fj : E + R(1 < j < s) are functions defined 


on a subset E + Ø of R, which are n times differentiable, and f® is the derivative of 
order p of the function f. Moreover, for every integer m’ Æ 0, the derivative of order 
p > 1 of the function f(x) = (x — A) , is given by fP (x) = m (m — 1)---(m' — p + 
p(m+p—1)! 


1)(@ — A)” -P, and when m’ = —m with m > 1, we get fP) (x) = (—1) / Di (x 
m—1)! 
àA) P. Now, applying the above formulas to the function H;(x), written under the 
1 S 
form H;(x) = —,; = II h;(x), with hj(x) = (x — àj)", we obtain 
I] (@-Aj)™ j=1j#i 
J=1, 74% 
Fi) (2) = 3 ji Ta] I ni?) (x), where êl] = {(p1, pos u Ps) ENP; pi +++ + 
ali j=1, j+i 


j i ;— 1)! 
Pi-1+ Disa +: + ps = k}. Since, hI (2) = (—1)?) (mj +p; — 1) (æ — Aj) Pi, for 


(mj; = 1)! 


every j (1 < j # i < s), we derive the following lemma. 


Lemma 2.8. For every k > land 1 < i < s, we have, 


OETA DA C p CH e-a o 


.— 1)! 
1... Ps Ah (m; — 1)! 


all 


Summarizing, through Expressions (16), (19) and (20), we can formulate the following 


result. 
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Theorem 2.9. Let A; (1 < i < s) be the roots of the polynomial P(z) = z" — aoz! —.--— Opa 
associated to the recursive relation (3), of multiplicities m4, m2,..., Ms, respectively. Suppose that 
for every root Ai the associated multiplicity m; > 2 (1 < i < s). Then, the analytic formula of the 
dynamic solution is given as follows, 


s mMi—l /mi—1 
r— < i— l Mi—a— — n 
S E O O wan, en 


for every n > r, where s(d, h) are the Stirling numbers of the first kind and H;(x) = 


and the H(z) are as in (20). 
Suppose that for every root A; the associated multiplicity m; = 1 (1 < i < r). Then, the 
analytic formula of the dynamic solution is given as follows, 


oS) a ; (22) 


for everyn >r. 


We illustrate Theorem 2.9 by the following special case. 
Special case: s = mı = M = 2. Let determine the dynamic solution o) when s = mı = 
mə = 2. By using Equation (21), we infer that, 


f 1 
s(d,h) Hox) n”? = Qi (n) + Q2(n), 
(Zren() One 


1 1 
1 = 
O(n) = D0 (> s(d,h) ( ) HË 2an) nl'X! = Biod} + BianX?, 


for i = 1, 2. Since 


=0 
s( = s(1,1) = 1 and s(1,0) = 0,a ee ne 
implies that 61, = H} 


o 
Oa), Bra = HOAT, bzo = HP (A2), Bor = HP (Arg) az?. 


Applying Equation (20), we obtain, 
—2 1 2 1 
Pio = Orda’ b11 = OSAA B29 = EO b21 = TESA 
Therefore, for s = mı = m2 = 2, the dynamic solution vÊ ) takes the form, 
—2 1 2 1 
(3) _ Ne4 Ne4 AR ++ d\n 
VU t n | , 
i = eu t Aea e eA Oaa 


for every n > 0. 


Let \; (1 < i < s) be the roots of the (characteristic) polynomial P(z) = z” — aọz"™t —- - - — 


ar—1, associated to the recursive relation (3), of multiplicities m1, M2,..., Ms, respectively. 
Without loss of generality, we set, 


Zı = {\;, root of P(z) with m; = 1}; Z = {A;, root of P(z) with m; > 2}. 
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Then, combining the two cases (21) and (22) of Theorem 2.9, we get the following general 
result. 


Theorem 2.10. Let A; (1 < i < s) be the roots of the polynomial P(z) = z maz OS 
associated to the recursive relation (3), of multiplicities m1, M2,..., Ms, P Then, the 
analytic formula of the dynamic solution is given as follows v’—) = 0) + YC), for every 
n > r, where 


mi—1 
1 7 
(r) n (r) _ (r) n? 
TUDE z waa- fr pai È Afin Jar 
i€Zı TI] Qi-An)™ i€Zq 
k=1, ki 
a mj — 1 1 
where Al") = S~ s(a, ny( , T with H,(«) = — ; 
d=h I] (@-A,)™ 
k=1, ki 
and the H(z) are as in (20). 
Since the analytic expression of of aa ) for p=1,. — 1, will be useful in the sequel, 
the result of Theorem 2.10 allows us to obtain, 
ge = ol”) d ae ws for every n > p, 
where 
O Mee () 1 (r) h| yn- 
Or r= >, 5 WP E py È Ra- ; 
iEZy II (A; = Ap) iE Zo h=0 
k=1, k#i 


m,-1 mi-1 h = mi—1 
> At? (n-p) = ` AY (>: (k) Mp) canbe written under the form b AG (n— 
h=0 


h=0 


mi—l /mi—1 
p)? = ye ( y =p é Ant LP 3 n". Therefore, we derive the following proposi- 


; k=0 \h=k 
tion. 


Proposition 2.11. Under the data of Theorem 2.10, for 1 < p < r — 1, we have, 


Ge au = >> Or + Se BD de, fi forn—p>r 


1€ Z1 1€ Zo 
( ) dP AP mi— l1 m—1 h ( ) 
r) ; í h—k r) a-k \ ok 
where bip = 5 and 3, = Glam! =i) 5 > (—1) (;) AnP ) Neg 
TI Qi —Ak)™ k=0 \ h=k 
k=1, kAi 
such that A?) = Y` s(d,h) (™ f ‘) HOED 0.) \-4 and H,(c) = ——— 
dah JI (@-Ax)™ 


where the HP) (x) are as in (20). 
116 


For illustrative purpose of Theorems 2.9 and 2.10, we consider the special useful two 
cases: s = 2 with mı > 2, mo > 2 and s = 3 with mı = məs = 1, mg > 2. For the first case 
we show that Z; = (and Z2 = {1, 2}. Thus, we have the following first corollary. 


Corollary 2.12. Special case s = 2: mı > 2 and mo > 2. Under the data of Theorem 2.10, we 
have, 


1) ATP mı—l /mı—1 m2—-1 /me2-1 
Un—p = (m = 1)! » ( >, qı(h, v) n* nt E aS )! Ji > l De 72 (h ») m3, 
" k=0 \ h=k k=0 \ h=k 


for everyn—p > r, with 1 < p < r — 1, where yı(h, k) = (=1)"-* (;) AS) ph-* ana 


k 3 
h r = : T TE -1 mı—a— — 
yalh, k) = (—1)*7 ay Ete = an(n a ame OCH) Pau 
d=h 
(r) ee mz —1 (mg—d-1) A 1 (k) 
Ash = Do s(d, h) d H5 (à2)A34, Hi(x) = — „and the Hy"? (x) 
sh I] -am 
k=1, ki 


are as in (20). 


For the second, we show that Z; = {1, 2} and Z = {3}. Thus, we have the following 
corollary, 


Corollary 2.13. Special case s = 3 with mı = mo = 1, m3 > 2. Under the data of Theorem 
2.10, we have, 


(r-1) ri” ‘i A3” i 


m | t Q3(n), 
Un—p (Ai = d2)(4 = A3 )™3 AY (A2 = A1) A2 = A3 )™3 A» a(n) 
—p m3—-1 /m3—1 
for n-— p > r, with1 < p < r — 2, where Q3(n) = oon se ( 5 y3(h, v) nfAR, 
3 | k=0 \ h=k 
h 3 . m3—1 ae eae 
ys(h, k) = (1)? AnA by s(a.n)(™ i Ja! RON and 
d=h 
H(x) = o such that the HP) (x) are as in (20). 
I] (@ — Ak) 


3. CONSTRUCTION OF THE INVERSE OF GENERALIZED VANDERMONDE 
MATRIX VIA THE ANALYTIC FORMULA (5) 


We have exhibited the close relation between the analytical form of sequences (3) and 
the generalized Vandermonde systems, through the equivalence of the two generalized 
Vandermonde systems (2) and (6). This section is devoted to the process of constructing 
the inverse . the generalized Vandermonde matrix, using the analytic formula (5) of the 
elements {uP Bae of the fundamental systems (7). 


We first introduce some useful notations, allowing us to introduce the generalized Van- 


dermonde matrix and to study its inverse. Let C-vector space of the polynomials C{z], 
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d 
and consider the derivation degree Dp(z) = r: Let po, pı, --, Pr—1 be in C[z] and 


f:C—C’, the valued vector function defined by, 


f(z) = (po(z), pı (z), eset ,Pr_1(z))", 
where (Y0, 71;---; 7-1)? denotes a vector column. We extend the derivation degree to 
the preceding vector function as follows, 
Df(z) = (Dpo(z), Dpr(z),--., Dpr-1(2))”- 

More generally, we have D® f(z) = (D™po(z), D™ pi(z),..., D™ pp_1(z))", for every 
k > 0, where D® = 1g, D® = D,..., D® = Do...oD, k times, for k > 2. 

Let A1,..., As be non-zero distinct s complex or real numbers, and s integers m1, ..., Ms, 
with r =m, +-+- + ms. Let C : C > C” be the valued vector function defined by, 

C(z) =(1,z,...,2" 1)", where z € C. 


For every k > 0, we consider the family of vector columns, 


ci = C(Ai) = (1, Ai, A2,- .., Pier k = 0 23) 
ch) = DMCA) = (0, i, 222, ..., (r — EMT? for k > 1. 


The generalized Vandermonde matrix of order r associated to the preceding family of vectors 
column (23), is given by, 


V= la, 0, E EN seo see), (24) 


According to [9], we have, 


S 


det V = (ù pea ( [O!L!... (m; — ni) Moare 
i=1 i=1 


j>i 
This expression shows that det V ¢ 0, because each \; # 0. Then, the generalized Van- 
dermonde matrix V has inverse V~!. 


Let built the process of the inversion of the generalized Vandermonde matrix (24), by 
utilizing the analytic formula (5) of the sequences of the fundamental system (7). For 
reason of clarity and conciseness, let consider the following useful notations of the two 
vectors column, 


T T 
= (Oromis bimi bs Or sm1) and A = (a0, Q1,- --,@r—1) ; 


where the p; j (1 < i < s, 0 < j < mi — 1) are the scalars of the analytic formula (5) and 
Qo, Q1, - - -, Ap—1 are the initial data given of sequence (3). Therefore, the linear system (5) 
can be written under the matrix equation, 


V-B=A. (25) 


Since the generalized Vandermonde matrix V is invertible, we derive, 


=Vit.a. 


For every p, n = 0,...,r — 1, we set A, = (0,...,0,1,0,...,0)" where 1 is located at the 
(p+ 1) —th position, or equivalently A, = (a0, a1,..., a1)", whith an = p,n: Then, the 
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(p + 1)-column of the inverse V~', of the generalized Vandermonde matrix (24) is given 
by, 


see gee mA (26) 
(-1) 


Furthermore, the scalar components of this column vector V,,,_;' are the solution of the 
linear system (5), defining the analytic formula of the linear recursive sequence (7), whose 
initial conditions are the entries of vector column A,. More precisely, for the sequence 
{uP Moy defined as in (7), and whose initial conditions are the entries of the column 
vector A,, the analytic formula is, 


ml 


= 2 As dnd | Ap. 


$1 j=0 


And by considering the vector column, 


Evie E oye Sd) 


we show that Expression (26) implies that the (p + 1) — th column, of the matrix V~! 
inverse of the generalized Vandermonde matrix, is given by, 


a: 
VS BO = (OBO ae Oi pO a 


1,0 


Therefore, the analytic formulas of the sequences {oP eer (0 < p < r — 1), defined as in 
(7), provides us the columns of the inverse of generalized Vandermonde matrix V given 
as in (24). In summary, we formulate the following fundamental result. 


Theorem 3.1. Let V be the generalized Vandermande matrix (24), related to the linear system 
(5), through the matrix Equation (25), namely, V -B = A, where B is the vector column B = 
(B10 ---, Blimi—1;-- -, Bs,0; <- -> Bs,ms—1)? obtained from the analytic expression (5) and A = 
(ag, Q1,---,Qp—1)/ is the vector related to the initial conditions of the sequence (3). Then, the 
inverse of the generalized Vandermande matrix is given by, 


v= BBO, B0], 
where ©) = (BP). eae mean OE. aN is the vector column Vo? = 
Vi. Ap with Ap = (a0,Q1,---,Qr-1)? with an = pn, namely, B®) = yT = 


VAA; 


Taking into account Theorem 2.2 and Theorem 3.1, as well as the results of Section 2, 
concerning the explicit analytic form the dynamic solution ol) we will give an explicit 
form of the vector column B®). Thus, we can establish an explicit form of V—! the inverse 


of the Vandermonde matrix V. 


In order to better grasp our process, let consider the special case of s = 2. Suppose that 
the characteristic polynomial of the sequence (3) is given by P(z) = (z — 1 )™ (z — 2)", 
for the sake of generality we take m1, m2 > 2. First, the related generalized Vander- 
monde matrix of order r = mı + mz, is given by, 


V= fa, Q, ios ny: C2, P, bsg cr), (27) 


where the c; (i = 1, 2), {P (1 < k < mı — 1) and Pa < k < mə — 1) are defind as in 


Expressions (23). 


Second, the analytic expression of each element {ul BN si where 0 < p < r — 1, of the 
associated fundamental system (7), can be written under the form, 


mı—l1 
P) — > ste) ee 
j=0 


Third, the associated column vectors 


mĐə—1 
| 5 se) AZ, for every n > 0. 
j=0 


(P) where 0 < p <r-—1, related to the former 


p) 


analytic formula of {vn }n>0, are given by, 


E. ey e 10s es 


In summary, we have the proposition. 


(28) 


T 
, B2 m2-1) "A 


Proposition 3.2. Under the preceding data, the inverse of the generalized Vandermonde matrix 
(27) is given by, 


v- = BO B®,..., 


(r=0] 


’ 


where the 


(P) are as (28). 


The previous fourth steps investigated to exemplify the preceding special case s = 2, 
allow us to formulate the following general algorithm for constructing the inverse of the 
generalized Vandermonde matrix, 


Step 1. The generalized Vandermonde matrix. Let \;,..., 
and the integer m1,..., Ms > land setr = mı +++- 
Vandermonde matrix of order r, is given by, 


às be in K (K = Ror ©), 
+ ms. The associated generalized 


EO cl), RAS clms—1)] 


(mi—1) 
Ci s. 9 


where the c; (1 < i < s) 


defined by (23). 

Step 2. Analytic form of the associated fundamental system. Let {oF AT 0<p< 

r — 1) be the fundamental system (7), whose characteristic polynomial is P(z) = [| (z — 
i=1 


-1 


Ai)" = 2" — aoz"™> — +++ — ar—1. Suppose that, for each p (0 < p < r — 1) the analytic 


É (Esp), 


Step 3. Vectors column _ y-t, For each p (0 < p < r — 1), we consider the vectors 
City Bey eae 


formula of {vl }n>0 is computed under the form yl 


column B®) = E E eee 


Step 4. The inverse of the generalized Vandermonde matrix V. The matrix V~! is given 
by V7? = BO, BY,..., BCD], 
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Special case: r = 3 and s = 2. Let \1, Az be in K and the integers mı = 1, mz = 2. Hence, 
we have r = mı + m2 = 3. In this case, we have cı = (1, A1, A7)", c2 = (1, A2, 43)? and 


ch) = (0, A2; 2A2)T. Therefore, the associated generalized Vandermonde matrix is, 


1 1 0 
V= |à à A2 
A2 A 2X2 
The characteristic polynomial associated to V is given by 
P(z) = (z — à )(z — à2)? = 2° — aoz? — az — a2. 
Let {uP ae, (0 < p < 2) be the fundamental system (7) related to the recursive relation 
Un = G0Un—1 + A1Un—2 + G2Un—3, for n > 3. Suppose that the analytic formula of each 


2 mi— l1 
{e jae, is given by, v®) = `> `> BO nd A?. Then, the inverse V~! of the matrix V 
i=1 \ j=0 
is as follows, 
0 1 2 
K By arg] 
V= 1650 830 Prol- (29) 
1 2 
La a 882] 


4. THE EXPLICIT ANALYTIC FORMULA FOR THE INVERSE OF 
THE GENERALIZED VANDERMONDE MATRIX 


4.1. First approach. The first approach of the dynamic solution of Subsection 2.2, allows 
us to exhibit explicit formulas for the entries of the inverse of the generalized Vander- 
monde matrix, using Propositions 2.5, 2.6 and Theorem 3.1. To reach our goal, we con- 
sider the following lemma, derived from Expression (8). 


Lemma 4.1. Let {vo }n>o (p = 0,...,r — 2) be the sequence of the fundamental system (7). 


s m,—-1 
Then, we have v?) = Sloe, BO nd ri, for every n > 0, where the che are written under 


i=1 
the form, 
pe = Ar—p— E ) + ar- w Fec E (30) 
such that the of ) are given as in (11). 
(p) _ (r= P a O EEE (r—1) 
Theorem 2.2 implies that vn” = ar-p—10p—1 + ar—-pUp—2 +: + ar—1Vp- p1; Moreover, 
s a 
Proposition 2.6 shows that EP = 5 y COn àz, for 1 < k < p+ 1, where the 
i=1 \ j=0 
oF are given as in (11). Therefore, we obtain, 
s mi—1 
vl) = 7 | Y (api? + pC +--+ aCe?) ni | ap, 
i=1 \ j=0 
for every p = 0,...,r — 2. Hence, the result follows. 
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By means of Theorem 3.1 and Lemma 4.1, we can exhibit the first result concerning the 
explicit expressions for the entries of the inverse of a generalized Vandermonde matrix. 


Theorem 4.2. Inverse of the generalized Vandermonde matrix YV. Let \1,...,s be non- 
zero distinct s real (or complex) numbers and s integers m1,...,ms. Let Y be the associated 
generalized Vandermonde matrix (24) of order r = mı + --- + Ms. Then, the inverse vo! of the 
matrix V, is given by, 


y- = pO, Bo 


’ ’ 


oid B070], 
such that B®) = (OU Ma iene E )T (0 < p < r — 1), namely, the 


’ 


S,Ms— 1 
explicit formulas for the entries by, (u,v = 1,...,r) of V-t = [buvli<u,v<r are as follows, 
eee for 1<u<m 

he meee for mt+1<u<m+m 
buv = i 

Bad ea eee wits, for my +My es + Ms—-1 + 1 sus my + M2 +: + Ms 
where the po are as in Expressions (9)-(10) and the an (p = 0,...,r—2) areas in Expressions 
(30), with the Cy are given by (11). 


For illustrating the result the former Theorem 4.2, we consider the following special 
case. 


Special case: r = 3 and s = 2. Let à, A2 be in two distinct complex numbers and the two 
integer mı = 1, mo = 2. Let V = [c1, c2, cP] be the associated generalized Vandermonde 
matrix of order r = mı + m2 = 3 defined by (23)-(24). Let consider the polynomial of 
degree 3 is given by, 
P(z) = (z= `) 2)? = 23 = (2A2 + à1)2? + (A2 + Dre) z= AA}. 

Thus, we have P(z) = 2° — aoz? — aiz — ag, with ag = 2A2 + Ay, ay = — (A2 + 21 A2) and 
az = 1X3. Let calculate the entries of V~' the inverse of the generalized Vandermonde 
matrix V. Taking into account s(0,0) = s(1,1) = 1 and s(1,0) = 0 and Equation (10) it 
follows 


’ 


1 
Ba = l(i A2), py = On, 22) e Ba = I AnA) 


: : (1] = -1 

And, using the formula (9) we derive the yọ (A1, A2) = w S yo (A1, A2) EE 

—1 
and P (Ai, A2) = yo: Therefore, we reach the formulas, 

1— Ag 

Gre, 1 @Q) = ne 
sei aa OEE Vi S ame OTST 

Equation (30) shows that By = RO e + aC À) and ao = ago), Since ay = —(A3 + 


212) and az = A, \3, then, a direct computation, applying the formula (11), allows us to 
establish, 

ty A+” 

Tt 2(A1 — 2)’ 


AB —2\1 A2 + A2? (0) _ 
Do = Az Bx ee aye bead = 
(A2 — A1) (Ai — 2)? dy — Ag 
Therefore, using the expression (29), we derive that the inverse V~! of V = [c1, ca, cP], is 
as follows, 


0 
p= 


| A2 —2)2 1 | 
(A2 = Ai)? A27) 2-A) 


y- = | —2\1À2 + `? 2X2 —1 
(A2571)? (An = Ar)? 2-A) 

| —Xj Ao + AY —1 
Maks aeea ae 


4.2. Second approach. Let consider the second approach of the dynamic solution of Sub- 
section 2.3, for establishing the inverse of the generalized Vandermonde matrix, using re- 
sults of Theorems 2.9, 2.10 and Theorem 3.1. To reach our goal, let consider the following 
theorem, derived from Expression (8). 

In order to identify the explicit formulas of the entries of the inverse matrix of a gener- 
alized Vandermonde matrix (24), let recall hereafter, the result of Theorem 2.9 concerning 
the dynamic solution (21), which can be adequate for its application in this subsection. 


Lemma 4.3. The expression of the dynamic solution vE” is given by the following formula, 


sS mi—1 
an ( Sbn” i every n > r, where 
i=1 h=0 


mi—l 


r—1) mi— 1 mi—d—1 = 
-iy Ser (h aeos o 


s(d, h) are the Stirling number of the first kind and Hie) (A;) are as in (20). 


Moreover, as we will make use of the formula (8), we have also the following property. 


s mi—l1 
Proposition 4.4. For d = 1,...,r — 1, we have op) = Ss" 5 Cni Xi, for every 
i=1 \ j=0 
n > r, where 


m1 
C9 = x4 > ybvate? (F) a4 (32) 
ae 
such that the ao are as in (31). 


In the aim to utilize Theorem 3.1 for constructing the inverse of the generalized Van- 
dermonde matrix (24), let now consider the following lemma, derived from Expression 


(8). 
Lemma 4.5. Let WP } aso (p = 0,...,r — 2) be the sequence of the fundamental system (7). 


s mi— l1 
Then, for every n > r, we have v®) = 5 e 5y ene ni j , such that ae are written under 


i=1 


the form, 


= = Ap—p— ice ) + ar- rele He tar E (33) 
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where the C are given as in (32). 


As a result of the above data, we get the following result regarding the explicit formula 
for the entries of the inverse of the generalized Vandermonde matrix. 


Theorem 4.6. Inverse of V. Let \1,...,As be non-zero distinct s complex or real numbers, 
and s integers m1,..., Ms. Let V be the associated generalized Vandermonde matrix of order 
r = mi +: + Ms, thus we have 


y- = RO, BO), C-D], 


gereg 


such that B®) = (3%, te er an pa, ae eA where the po are given Ex- 
pressions (31), with HCD (Ai) are as in (20), and the pe (p = 0,...,r — 2) are given by 


j 


Expressions (33), with the a are as in (32). More precisely, the explicit formulas of the entries 


buy (u,v =1,...,7r) of V-t = [bwwli<uv<r are as follows, 


By for 1l<u<m 


. Pies a a for m+1<u<mi+my 


E ap JOY Mmi + m +e + Ms- +1 <u < m +m tt ms 


s,u— (mi +m +- +m 


where the Bee are given Expressions (31), with HERG are as in (20) and pe (p = 


(d) 


0,...,7 — 2) are given by Expressions (33), with the OF are as in (32). 


In order to better understand the general case of Theorem 4.6, we present below, a 
significant special case. 


Special case: r = 4 and s mı mə 2. Let Ai, Az be in two distinct complex 
numbers and the two integer mı = 2, m = 2. Let V = [c1, cf) C2, cP] be the associated 


generalized Vandermonde matrix of order r = mı + m2 = 4 defined by (23)-(24), namely, 
1 0 1 0 

Ar Ar Ag Ag 

Ne DAS oe DA? 

A3 3A A3 BNE 

Consider the polynomial of degree 4 given by, 


v= 


P(z) = (z—1)?(z— 2)? = 24 — aoz? — az? — azz — ag 

where ag = 2(A1 +2), ay = — (A2 +41 A2 + AP), ag = 2(A2A1 + A? A2) and a3 = —(A?3). 
Using result of the special case illustrating Theorem 2.9, we have, 

a(S) 2 PMA OY AUN BG) PAIN, S E 

BOT Ar = ABALA PE (Ar = AAA O (Ar = ABALA a (Ap = A2)8ALAQ” 
From Equation (32) and Equation (33), and because a, = —(A3+4A1A2+A7), a2 = 2(A3A1+ 
A22) and a3 = —(A7\3), a direct computation shows that, 
323 — ALAS _4(0) —A2A3 + ALAS (0) —3A3A2 + ATAg. 

1,1 


(Ar — à2)8A1A2 P (Ar — A2)8 Ad A279 29 (A, = A2)8ALA2’ 
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A(0 
BY) = 


o _ ZATA +AA gy GATAQ aa) _ ALAR + QATAR A3, 
tT Or = A2)8A1A2 PP OA = Ag)BALAg ht (Ay — A2)3A1A2 
gy 7 6A2A2 a _ ABA — 2AZAZ + AT. a2 2 3A1A3 + 3BATAQ | 
2,0 (Ài = A2)5A1A2 2,1 (ài _ A2)3A1A2 re LO (Ai = à2)5A1A2 
32) =, —r,r3 + 2A3 — AZAZ. 5,(2) = —3)\,A3 = 3A? r2 . 32) = dr? do + AA3 — 2A 
at (A1 = ABAAA 8 Ad = AAA 7%? (A1 — A2)3A1A2 
Therefore, using Theorem 3.1, we derive that V~! the inverse of V = [c1, cf) c2,€ 
8 1 
as follows, 


8s 


(1)) is 


SAGAS = Air} —6A2 03 3A1\3 + 3A2A2 —2\1 A2 
T ENS H AAF AA H 2A2AZ = AG AMAZ POS ae AA A2 
~ K |383 + APA 622 —3A,A2 — 3422 2AA |” 


—AfA2 + AZAR ABAQ — QAPAB HAT — ARA + AAZ — 2A A? = AA 
where K = (` — A2)?A1A2. 


5. ANALYZE AND DISCUSSION 


In view of their numerous applications, the usual Vandermonde systems and gen- 
eralized Vandermonde systems, have been largely studied in the literature (see for ex- 
ample [6,7, 10, 15, 17,20] and references therein). In order to resolve them, the topic of 
the inverse of associated Vandermonde and generalized Vandermonde matrices contin- 
ues to attract a lot of attention, and has been the subject of numerous research papers. 
Therefore, various approaches have been implanted for succeeding this inversion (see 
for example [2,7, 11, 15, 16, 18,20] and references therein). Especially, the technique of 
LU factorization of matrix (see [12, 15,19]), has been examined in [12], for inverting the 
generalized Vandermonde matrices. 


In our study, we have exploited the fact the resolution of the Vandermonde and gen- 
eralized Vandermonde systems (1)-(2), is linked to the determination of the analytical 
formula (5) of the linear recursive sequences (3). Recently, such important relation has 
been exploited in [2] to develop a method for inverting the Vandermonde and general- 
ized Vandermonde matrices. More precisely, Expression (13) from [2, Theorem 2.9] has 
been utilized in the aim to establish some explicit formulas of the 8; ; (see [2, Proposition 
4.1]). 


Regarding our results, we approached the inversion of the generalized Vandermonde 
matrix via the computation of the scalars 6; ; or ĝi; (1 < i < s, 0 < j < mi — 1), by 
considering other explicit formulas for analytical expressions of the linear recursive se- 
quences (3), and through another approach based on properties of linear algebra. More 
precisely, such explicit analytical expressions are applied to the fundamental system (7), 
which makes it possible to construct the vectors columns of the inverse of the generalized 
Vandermonde matrix V~!. This represents a new procedure, in the algorithmic construc- 
tion of the matrix V~', using Theorem 3.1. It is important to specify that, when all the 
roots of the polynomial P(z) = 2” — agz’~! — --- — ar—ı are simple Expression (16) and 
(22), permit to exhibit the inverse of generalized usual Vandermonde matrix as shown 


in [2, Theorem 3.3]. And by using our new process this result can be recovered easily. 
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6. CONCLUSION AND PERSPECTIVES 


This study presents results regarding some explicit formulas for the entries of the in- 
verse V`! of the generalized Vandermonde matrices. Our process based on the explicit 
analytic formula of fundamental system (7), and an algorithmic technique, for construct- 
ing the columns of V-t, In the best of our knowledge, the formulas of the entries of vol 
are not current in the literature. 


We do not claim that this study is a complete presentation of all methods for approach- 
ing the the inverse V~! of the generalized Vandermonde matrices. However, the material 
presented here is a reflection of our approaches concerning the inverse V~! of the gen- 
eralized Vandermonde matrice, through the properties of fundamental system related to 
the associated to difference equations (3). 


Finally, as there are other explicit analytic forms (5) of the general term of the linear 
recursive sequence (3), this opens up possibilities for obtaining other explicit formulas of 
the entries of the inverse V~! of the generalized Vandermonde matrices. 
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